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We estimate the energy loss distribution and investigate the quenching of hadron spectra in 
] ultrarelativistic heavy-ion colhsions due to the collisional energy loss of energetic partons from hard 

. parton collisions in the initial stage. 
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CO ', In the initial stage of ultrarelativistic heavy-ion collisions energetic partons are produced from hard collisions 
• between the partons of the nuclei. Receiving a large transverse momentum, these partons will propagate through the 
fireball which might consist of a quark-gluon phase for a transitional period of a few fm/c. These high-energy partons 
will manifest themselves as jets leaving the fireball. Owing to the interaction of the hard partons with the fireball 
qJq these partons will lose energy. Hence jet quenching will result. The amount of quenching might depend on the state 
\^ , of matter of the fireball, i.e., quark-gluon plasma (QGP) or a hot hadron gas, respectively. Therefore jet quenching 
has been proposed as a possible signature for the QGP formation [1] . 

In order to use jet quenching as a signature for the quark-gluon plasma, the energy loss of hard partons in the QGP 
has to be determined. First the energy loss due to collisional scattering was estimated [2-6]. Using the Hard-Thermal- 
\ Loop (HTL) resummed perturbative QCD at finite temperature [7] , the collisional energy loss of a heavy quark could 
be derived in a systematic way [8-11]. From these results also an estimate for the collisional energy loss of energetic 
' gluons and light quarks could be derived [12], which was rederived later using the Leontovich relation [13,14]. Later 
also the energy loss due to multiple gluon radiation (bremsstrahlung) was estimated and shown to be the dominant 
O ^ process. For a review on the radiative energy loss see Ref. [15]. Recently, it has also been shown [16] that for a 
D , moderate value of the parton energy there is a net reduction in the parton energy loss induced by multiple scattering 
i-S^ ' due to a partial cancellation between stimulated emission and thermal absorption. This can cause a reduction of the 
^ [ quenching factor due to radiative processes. 

Unfortunately, jets, requiring very large initial parton energies, are rare events, which are difficult to observe. 
However, quenching of hard partons will also affect hadron spectra at high transverse momenta pj^. Indeed, first 
results from RHIC have indicated a suppression of high-p^ spectra [17]. The consequences of jet quenching on hadron 
spectra have been calculated, for example, in Refs. [18,19]. Here only the radiative energy loss has been taken into 
account. The purpose of the present paper is to estimate the quenching of hadron spectra due to the collisional energy 
loss of partons in the QGP. As we will see, this contribution can be of the same order as the radiative quenching. 

We will follow the investigations by Baier et al. [18] and Miiller [19], using the collisional instead of the radiative 
parton energy loss. Following Ref. [18] the p± distribution is given by the convolution of the transverse momen- 
tum distribution in elementary hadron-hadron collisions, evaluated at a shifted value p± + e, with the probability 
distribution, -D(e), in the energy e, lost by the partons to the medium by collisions, as 
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Here Q{p±) is suppression factor due to the medium and the total energy loss by partons in the medium is 



A£; = y eD{e)de. (2) 



We need to calculate the probability distribution, D(e), that a parton loses the energy, e, due to the elastic collisions 
in the medium. This requires the evolution of the energy distribution of a particle undergoing Brownian motion. The 
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operative equation for the Brownian motion of a test particle can be obtained from the Boltzmann equation, whose 
covariant form can be written as 

jfd^D{x,p)=C{D}, (3) 

where p is the momentum of the test particle, C{D} is the collision term and D{x,p) is the distribution due to the 
motion of the particle. If we assume a uniform plasma, the Boltzmann equation becomes 



dD _ C{D} _ /aD\ 
'dt ~ E ~ \ dt) 



(4) 

coll 



Wc intend to consider only the elastic collisions of the test parton with other partons in the plasma. The rate of 
collisions w{p, k) is given by 

w{p,k)= ^ w^{p,k), (5) 

where represents the collision rate of a test parton i with other partons, j, in the plasma. The expression for 
can be written as 

w\p,k)=-ijj J^fj{q)Vreia^ , (6) 

where jj is the degeneracy factor, v^ci is the relative velocity between the test particle and other participating partons 
j from the background, and is the associated cross section. Due to this scattering the momentum of the test 
particle changes from p to p — k. Then the collision term on the right-hand side of (4) can be written as 

= J d^k[w{p + k,k)D{p + k)-w{p,k)D{p)] . (7) 

where the collision term has two contributions. The first one is the gain term where the transition rate w{p + fc, k) 
represents the rate that a particle with momentum p + k loses momentum k due to the reaction with the medium. 
The second term represents the loss due to the scattering of a particle with momentum p. 

Now under the Landau approximation, i.e., most of the quark and gluon scattering is soft which implies that the 
function w{p,p') is sharply peaked at p « p', one can expand the first term on the right-hand side of (7) by a Taylor 
series as 



d 1 

w{p + k,k)D{p + k) w{p,k)D{p)+k - —{wD) + -hkj {wD) + (8) 

C/p Z OPiOPj 

Combining (4), (7) and (8), the Fokker-Planck equation 
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k . A + ifc A; ^ 
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{wD) (9) 



is obtained, which describes the equation of motion for the distribution function of fluctuating macroscopic variables. 
We consider, for simplicity, the one dimensional problem, for which (9) can be written as 

Here the transport coeflacients are given as 

T,{p)= J d^kw{p,k)k=^^ = {F), 

T2{p) = Skw{p,k)k^ = . (11) 

Now the work done by the drag force, 71, acting on a test particle is 
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-dE = (F) ■ dL = Ti{p) ■ dL , (12) 
which can be related to the energy loss [9,12] of a particle as 

dE 

--^=Ti{j>)^pA, (13) 

where A is the drag coefficient, which is almost independent of momentum p [3,20]. The drag coefficient is a very 

important quantity containing the dynamics of clastic collisions. In principle, A{p,f^, may be determined from the 
kinetic theory formulation of QCD through the application of the fluctuation-dissipation theorem [21], but that is 
indeed a difficult problem. As discussed in Refs. [3,20,22], the drag coefficient is expected to be largely determined 
by the properties of the bath and not so much by the nature of the test particle. Then one can use the average of the 
drag coefficient given as 

{A{p,t)) = A{t) = ^-i . (14) 

The quantity T2 can be related to the diffusion coefficient in the following way: 

T2{p) = = pAp « AT^ = Vf , (15) 

where we have approximated p by the temperature T of the system and the drag by using the Einstein relation, 
TiT « Dp, assuming that the coupling between the Brownian particle and the bath is weak [21]. 
Combining (10), (13) and (15), we find 

-=A-ipD) + Vp—, (16) 

which describes the evolution of the momentum distribution of a test particle undergoing Brownian motion. 
Next we proceed with solving the above equation with the boundary condition 

D{p,t) S{p-po) . (17) 

The solution of (16) can be found by making a Fourier transform of D{p, t), 

Dip, t) = ^ D{x, t) e'P-dx , (18) 

where the inverse transform is 

/+00 
D{p,t) e-'P^'dp . (19) 
-00 

Under the Fourier transform the corresponding initial condition follows from (17) and (19) as 

D{xo,t = to) = 6-'^°="° (20) 

where x = xq &t t = to is assumed. 

Replacing p i-^ and ^ ix, the Fourier transform of (16) becomes 

'4+Ax^ = ^Vpx^D. (21) 
ot ax 

This is a first order partial differential equation which may be solved by the method of characteristics [23]. The 
characteristic equation corresponding to (21) reads 

dt ^ ^ db 

1 Ax Vpx^b ■ ^ ' 

Along with the boundary condition in (20) and A{to) = Vpito) = 0, the solution of (21) can be obtained as 
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D[x, t) = exp 
Combining (18) and (23) yields 



exp 



—ipo x{t) e 



- f A(i')dt' 



exp 



(23) 



(24) 



It is convenient to integrate over xq instead of x and to substitute the solution of the first pair of (22) into the above 
equation, leading to 



D{p,t) = 



exp ( /* A{t') dt' 
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/ + 00 
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(25) 



Using the standard form 
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one can easily perform the integration in (25), resulting in 
D{p,t) = ^ 



-{z-K,f/K2 



(26) 



TT (4 /*I?F(t') exp 2fA{t")dt" dt'^ exp (^-2 f A{t") dt"^ 



X exp 



(a f Vpit') exp 2f'A{t")dt" dt'^ exp (^-2 f A{t") dt"^ 



(27) 



For relativistic particles, p = E, (27) can be written as 



D{E,L) = 



exp 



W{L) 



(28) 



where W(L) is given by 



W{L)= (^aJ^ VF{t')exp 2 j' A{t")dt" 



dt' 



exp 



! / A{t')dt' 
Jo J 



(29) 



The energy loss of partons in the QGP due to elastic collisions was estimated in Ref. [12] and we used the expression 
averaged over parton species as 



2"//2(6+"/) Q_92 . 



^ET 



rrio 



(30) 



where n/ is the number of quark flavours, ag is the strong coupling constant, nig = + Uf /6)gT/3 is the thermal 
gluon mass, and E is the energy of the partons. Following (14), we can now estimate A for different quark flavours and 

gluons at the energies (temperatures) of interest. For averaging over the momentum the Boltzmann distribution was 
used. The time dependence of the drag coefficient comes from assuming a temperature T{t) decreasing with time as 
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the system expands, according to the Bjorken scahng law [24] T{t) ~ tQ^^To/t^^"^, where to and Tq are, respectively, the 
initial time and temperature at which the background of the partonic system has attained local kinetic equilibrium. 
Since the plasma expands with the passage of time, we have used the length of the plasma, L as the maximum time 
limit for the relativistic case (w ~ 1). In Fig. 1 the drag coefficient in the QGP phase of the expanding fireball is 
shown as a function of time, where we have chosen the parameters Tq = 0.5 GeV, to = 0.3 fm, and as — 0.3. 
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Figure 1: The drag coefficient A{t) in an expanding QGP. 



In Fig. 2 we show the probability distribution, D{E, L) given in (28), as a function of energy E, choosing Eq — 1 
GeV. The peak of the probability distribution is shifted with passage of time (or distance travelled) indicating the 
most probable energy loss due to elastic collisions in the medium. 
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Figure 2: The energy-loss probability distribution, D{E,L) as a function of energy E for a given distance L. 

After traversing a distance L, the mean energy of the parton due to the elastic collisions in the medium can be 
estimated as 



'E) = 



ED{E,L)dE ^ Eoe 



- A(t')dt' - A(t')dt' 

Jo ' ' — rn± e •'0 



(31) 



where Eq = m±^ = y/p^^ + ra^ at the central rapidity region, y = 0. The average energy loss due to elastic collisions 
in the medium is given by 



AS 



, - A(t') dt' 

1 — e Jo 



For the massless case, the average Ap^ can be written as 

A ( ^ - A(t')dt' 



(32) 



(33) 
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Figure 3: The effective shift of the scaled transverse momentum Ap±/p± as a function of distance L. 



In Fig. 3 we display the scaled effective energy loss of fast partons due to elastic collisions in the medium as function 
of the distance L. It is found that the effective energy loss, defined as the shift of the momentum spectrum of fast 
partons, depends on the transverse momentum p± and the distance L. For given p± the effective energy loss is 10% 
after traversing a distance of Ifm and around 50% after 10 fm. The effective energy-loss scales linearly with p± for a 
given L. This clearly reflects a random walk in p± as a fast parton moves in the medium [19] with some interactions 
resulting in an energy gain and others in a loss of energy. 

We assume that the geometry is described by a cylinder of radius R, as in the Boost invariant Bjorken model [24] 
of nuclear collisions, and the parton moves in the transverse plane in the local rest frame. Then a parton created at 
a point r with an angle (j) in the transverse direction will travel a distance [19] 



2 • 2 

r sm ( 



Nl/2 



(34) 



where cos0 = v • f ; v is the velocity of the parton and r = |f|. 

The parameterization of the p± distribution [19,25] which describes the first RHIC hadroproduction data for mod- 
erately large values of p± has the form 



d^'p^L \ Pa 



(35) 



where « 8 and pq — 1.75 GeV. The quenched spectrum convoluted with the transverse geometry of the partonic 
system can be written from (1) as 



dN 
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(36) 
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Figure 4: The quenching factor Q{p±) as a function of transverse momentum p±. 



In Fig. 4 our numerical results of the quenching factor Q{p±) due to elastic collisions in the medium and its 
dependence on p± are shown. It is found that the quenching becomes stronger at lower p± and gradually decreases 
with higher p± in agreement with experiments [17]. As discussed above the quenching factor obtained here is quite 
similar to a random walk as well as the Bethe-Heitler scaling law given phenomenologically in Ref. [19]. 

Comparing the quenching factor due to elastic collisions, shown in Fig. 3, with the one coming from the radiative 
energy loss [19] we observe that both are of similar magnitude. The radiative quenching factor varies for p± from 6 
to 16 GeV between 0.17 and 0.24 depending on the model assumptions, whereas the collisional quenching factor lies 
between 0.22 and 0.26 for the same momentum range and for our choice of the parameters, i.e., Tq — 0.5 GeV, to — 0-3 
fm, and as = 0.3. Apart from uncertainties in these parameters let us also have a look at some of the assumptions 
made in this work which may modify the quenching factor. First, as discussed above, the momentum dependence of 
the drag coefHcent, containing the dynamics of the elastic collisions, has been averaged out. A major advantage of 
this is the simplicity of the resulting differential equation. Of course, this simplification can lead to some amount of 
uncertainty in the quenching factor. Secondly, the entire discussion is based on the one dimensional Fokker-Planck 
equation and the Bjorken model of the nuclear collision, which may not be a very realistic description here but can 
provide a very intuitive picture of the problem. However, extension to three dimension is indeed an ambitious goal, 
which may cause that many of the considerations of the present work will have to be revised. Within the limitation of 
our simplified model the contribution of the collisional energy loss of partons in the QGP to the quenching of hadron 
spectra in ultrarelativistic heavy-ion collisions is found to be very important and cannot be neglected. In particular, if 
the radiative energy loss is suppressed at intermediate energies, the collisional quenching would explain the observed 
decreasing of the quenching factor with increasing transverse momentum. 
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